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ABSTRACT 

In this paper, we develop a new gravitational lensing inversion technique. While traditional approaches 
^^ ' assume that the lensing field varies little across a galaxy image, we note that this variation in the field 

^^ , can give rise to a "Flexion" or bending of a galaxy image, which may then be used to detect a lensing 

^3 ' signal with increased signal to noise. Since the significance of the Flexion signal increases on small 

^N I scales, this is ideally suited to galaxy-galaxy lensing. We develop an inversion technique based on the 

"Shapelets" formalism of Refregier (2003). We then demonstrate the proof of this concept by measuring 
a Flexion signal in the Deep Lens Survey. Assuming an intrinsically isothermal distribution, we find 
from the Flexion signal alone a velocity width of Vc = 221 ± I2km/s for lens galaxies of r < 21.5, subject 
to uncertainties in the intrinsic Flexion distribution. 
Subject headings: galaxies: halos — galaxies: structure — gravitational lensing 
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^ . 1. INTRODUCTION 

, 1^1. . The past several years have seen an explosion in the analysis of weakly gravitationally lensed images of galaxies by 
■^ ' galaxies (e.g. Brainerd, Blandford & Small 1996; Hoekstra, Yee, & Gladders, 2004), clusters (e.g. Small et al. 1997; 
^ [ Wittman, 2001; Gray et al. 2002; Taylor et al. 2004), and Large Scale Structure (e.g. Wittman et al. 2002; Hoekstra et 
\^ . al. 2002; Jarvis et al. 2003; Brown et al. 2003; Pen et al. 2003; Bacon et al. 2003) producing an unprecedented glimpse 
f^ ■ into the underlying matter distribution of the universe. Galaxies-galaxy lensing, in particular, has proven a fertile testbed 
^^ , of our understanding of structure formation. 

C3 ■ Dynamical estimates of galaxy masses are subject to uncertainties about whether a system is relaxed, and are lim- 
ited by the lack of luminous dynamical probes beyond a few tens of kpc from the center of a galaxy. Gravitational 
Q^, lensing techniques, on the other hand, provide an accurate way of computing the surface density of distant objects 
I ■ without recourse to dynamical estimates, since the distortion of images is dependent only upon the potential field of 
O ' a lens and not upon either its composition or dynamics. The success of these endeavors has been such that there are 
a number of ongoing surveys of weak lensing fields (e.g. Wittman et al. 2002, hereafter DLS; CFHT Legacy Survey, 
^ http://www.cfht.hawaii.edu/Science/CFHTLS-DATA/) and instruments (e.g. Advanced Camera for Survey, Clampin 
et al. 2000; Dark Matter Telescope; Supernova/ Acceleration Probe, Rhodes et al. 2004) largely designed around the 
acquisition of large high-quality lensing data. Moreover, there remains a potential bounty of information to be found in 
existing and ongoing survey data (e.g. SDSS, York et al. 2000; Great Observatories Origins Deep Survey, Giavalisco et 
al. 2004; Medium Deep Survey, Ratnatunga, Griffiths & Ostrander 1999). Given the difficulty and expense in collecting 
C^ ■ high-quality surveys of lensed galaxies it is imperative that we extract as much information as possible from them. 

One of the great advantages of studying weak lensing fields is that the physics underlying gravitational lensing is well 
understood (see e.g. Blandford & Narayan, 1992; Kaiser & Squires, 1993; Kaiser, Squires & Broadhurst, 1995, hereafter 
KSB; Mellier 1999; Bartelmann & Schneider, 2001 for reviews of weak lensing). The KSB approach provides the standard 
for most current analyses of weak lensing, and is based upon estimating the ellipticity of a lensed galaxy as a probe of 
the local shear field. By measuring the ensemble properties of ellipticity and orientation for a number of sources, one can 
make a determination of the properties of the lens. 

While the currently applied approaches have done an excellent job in estimating the matter distribution of gravitational 
lenses, they ultimately measure only the ellipticity of an image (the second moments), and thus potentially drop significant 
information from sources with significant substructure. In order to improve on this, we plan to extend the work of Goldberg 
& Natarajan (2002) who suggested that second order effects in gravitational lensing fields may give rise to an octopole 
moment in the light distribution, which expresses itself as "flexion" in the image. 

The method of Goldberg & Natarajan ultimately relied on a very complicated form of the 2nd order shear operator, 
which made a practical inversion difficult. Instead, here we will cast this approach into a "Shapelets" formalism (Refregier, 
2003; Refregier & Bacon, 2003), a novel approach to both image and lensing analysis. Rather than analyze image 
shapes according to their multipole moments, shapelets methodology decomposes images into combinations of Hermite 
polynomials. As a reminder to the reader, the reduced Hermite polynomials are the eigenfunctions of the simple harmonic 
oscillator in quantum mechanics, Bn{x). These functions have a number of useful properties, including orthogonality and 
a Gaussian factor, which localizes the function. 

In § 2, we begin by reviewing basic properties of the "forward" problem in weak gravitational lensing, and define our 
notation within the present work. We then proceed to introduce the second order term in the lensing operator, and then 
discuss a particularly useful approach to the second order problem using the shapelets formalism. In § 3 we discuss a 
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2 Second Order Lensing 

practical technique for inversion of the density gradient signal. In § 4 we provide a proof of concept by measuring the 
second order shear signal from Galaxy-Galaxy lensing in several DLS shear fields. We conclude with a discussion of future 
prospects. 

2. SECOND ORDER LENSING 

2.1. Review of Weak Lensing Formalism 

We begin with a brief review of weak lensing theory. An excellent discussion of this material can be found in Bartelmann 
& Schneider (2001), from which we borrow our conventions. Imagine that we are observing a set of extended sources at 
an angular diameter distance, Ds which are gravitationally lensed by a mass distribution at an angular diameter distance, 
Di. We thus define a dimensionless surface density, the convergence, k, such that: 

»(x) . £^i^ . (1) 

where E is the surface density of the lens. Dig is the distance between lens and source and x represents the image 
coordinates as seen by the observer (neglecting a constant coordinate transform). 

The convergence may be thought of as a source term for a potential, ip{(^)i ^^id related via a Poisson-like equation: 

V2^(x) - 2«:(x) , (2) 

where all gradients and divergences are calculated in two dimensions. 

Since lensing conserves surface brightness, a mapping from foreground to background coordinates is sufficient to deter- 
mine a background brightness map from a foreground one (or vice-versa) provided a full knowledge of the geometry of 
the system (cosmology plus the redshifts of the source and lens) and mass distribution of the lens. Thus, we may expand 
around the origin to determine a deprojection operator on a foreground light distribution, which yields the amplification 
matrix, 

^ '- 9x ~ r^ dx^dx,) -\ -12 l-« + 7i j ■ ^^' 

where x' are the coordinates that the lensed image would have in the absence of lensing. Note that for convenience we set 
the origins of both the foreground and background coordinate systems to be the centers of light in their respective planes. 

Equation (3) is the first term in a Taylor series expansion of the distortion operator. The term 7 is a complex shear 
term, representing the anisotropic part of the distortion, with 7 = 1716^*"^, and the real and imaginary parts being denoted 
with the subscripts, "1" and "2" respectively, as per convention. 

If we assume that the convergence and shear field were constant along the scale of the lensed image, the brightness map 
of a galaxy would transform with the simple relation: 

x[ = A.jXj . (4) 

2.2. Expansion to 2nd Order 
Of course, the field is not a constant. Thus, we may imagine expanding the field to second order such that: 

Xj^ — ^ij'^j 1 "^J-^ijk-^j^k \^ ) 

where 

As Kaiser (1995) has pointed out, there are a number of relationships in the derivatives of k and 7 which greatly simplify 
the construction of D: 

f Ks \ _ f 71,1 +72,2 

Thus, 



K,2 J \ 72,1 -71,2 ' ^ ^ 



A,. ^ I -^^i-^- :5:' I (8) 



Dij2 - 



-72,1 —72,2 
-72,2 271^2-72,1 



Now, let us suppose that there exists some background brightness field, /(x'), such as the brightness of some background 
galaxy. Since lensing preserves surface brightness, we may write the apparent foreground brightness field as: 

/(x) = /'[x'(x)] . (9) 



We can write the brightness as: 
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/(x)~/'(Ax+-Dx<»x) 

= /'[x+(A-i)x+-Dx(^x] 



(10) 



Or, expanding as a Taylor series: 



/(x) 



[-^ ^ /ij-^j ~r c\ ^jk-^j-^k 
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dxi 



[A - I)ikXk + -DikiXkXi 



(11) 



(^ ^ I)jkXk + -DjklXkXi 



d^ 



dxidx-i 



/'(x) 



Expanding out the terms on the first line the expressions are first order in 7, and (up to) 2nd order in position, while the 
terms on the second line are (at least) 2nd order in position, and 2nd order in 7. If the scale on which 7 varies is typically 
smaller than the image size, the terms on the second line will necessarily be smaller than the 2nd order contribution on 
the first, and thus we will consider only the first line. 

Thus, in the limit of a smoothly varying, weak shear field, we have the relation: 



/(x) 



\-^ ^ )ij-^j ~r cy^ijk^j^k 
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dxi 



/'(x) . 



(12) 



2.3. 



Image Analysis Using Shapelets 

While equation (12) defines a linear transformation on an image due to lensing, in practice it is fairly complicated 
to actually perform such a transformation (or the inversion) in generality. However, we may simplify this problem 
considerably by decomposing the image into basis coefficients. 

Following Refregier & Bacon (2003), we expand the image into Reduced Hermitc polynomials, or "Shapelets." The 
light distribution, /(x, y) of a galaxy is expanded as a combination of two-dimensional basis functions: 



/(x) 



Y^fn^B^m 



(x) 



(13) 



where 

and where /3 is a scaling factor, and 

such that: 



are the reduced Hermite polynomials: 

0„(x) = [2"7ri/2„!J n„{x)e~'^ 



W' - 2xH' + 2nHn = . 



(14) 



(15) 



(16) 



Since the lowest order polynomials resemble Gaussian light profiles, Refregier (2003) demonstrates that for typical HST 
images, convergence can rapidly be reached using a few x 10 coefficients. 

It should be noted that though this expansion is not the same as a multipole expansion, it has many similar properties. 
For example, an image with /20 — a^, and /( 



02 



b^ < c? will look like an ellipse with a Gaussian radial profile. Likewise, 



any combinations for which n + m is even produces an image which is symmetric with respect to 180 degree rotations 
(see, e.g. Figure 2 in Refregier, 2003). The canonical pictures of both spiral and elliptical galaxies have precisely this 
symmetry, and thus, it is possible that typical galaxies can be reconstructed almost exclusively from even moments. 
If, taking our example from the SHO in quantum mechanics, we define: 
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and expand Equation 12 to first order in 7, then we get: 



./(x) 



1 + (A - I)ijX.jPi + -UijkXjXkPt 
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/'(x) 



(17) 
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Again following Refregier & Bacon (2003), we expand our light function in Shapelets and can apply the distortion 
operator as a combination of raising and lowering operators. Note that we have implicitly defined our coordinates such 
that the (3 parameter used in Refregier (2003) is equal to 1. Recall, 

ii = -7={ai + a\) Pi ^ -^{a, - a]) (19) 

and can operate on a Shapelet (e.g. with i = 1) as: 

ai |0n m) = \Al|0n-l m) (20) 

a\ \cj)n m) =\/n+l \(t>n+l m) (21) 

and similarly for the y-axis. 

Now, the first thing to note is that with the addition of the D operator, there are odd combinations of x and p, meaning 
that there can be coupling between An + A?n = odd modes. Contrary to ordinary weak lensing, in which an image is 
lensed symmetrically around its center of light (and thus, the centroid remains fixed), second order effects will cause a 
shift in the center of light. We must thus compute the shift in the center of light. 

Because of the inherent symmetries in the linear lensing operator, the only contribution to the centroid shift comes 
from the second order operator: 

{xi) = -Dijk j d^x xiXkXk— (22) 

Expanding this out and integrating, we find that the center of light in the foreground will be shifted by: 

(a;i) = -(a;?) {\diii + D2i^ - {x[x'^){2Dn2 + D222) - {x'i) ^^^122) (23) 

(3:2) = -(a;?) (^^211 j - (a;'i4> (2^212 + i^iii) - (4") [^^222 + Dn2 

where the second moments of the light distribution can be computed from the foreground field (since changes in the 
moments from background to foreground will necessarily be 2nd order in 7, and thus negligible). 
A translational shift can be given by the operator: 

f,^^{al-a,) (24) 

Combining all of this yields the following lensing operator: 

/(x) ^ (1 + «X + 7,5f ) + 7.,,.5|f )/(x') (25) 

where K is the convergence operator: 

K = l + ^{a\^ + al''-al-al) (26) 

and Sl is the linear shear operator defined by Refregier (2003): 

S['^ = l{ar~al'-&1 + &1) (27) 

S-W ^ a\al - aiaa . (28) 

It is straightforward, albeit tedious, to derive the explicit form of the second order operator, S^"^' from equations (18) and 
(22-24). We write out the explicit 2nd order transformation as: 



:.(2) 



s3 I ;i. / /I _ o /Cr I 1 o ^ -r,^ \ I _L a / -„,, \ A a ^ ^,, \ ;;t i ^t 



S\^' = —^ -2a( + ai [4 - 2N + 12 < XX >]+ 8 < xy > 0.2 ~ 8 < xy > al + al [6 + 2N - 12 < XX >] (29) 
+2af 



c(2) _ 1 

O-in 



'' "4V2 
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(31) 

(32) 



where NBnrn{x,y) = nBnm{x,y), and similarly for M. 

Interestingly, all of the terms in the second order shear couples An + Am =odd coefficients, while the first term shear 
and convergence operators couple An + Am = 2. In other words, we may consider only the "odd" (meaning n+m=odd) 
moments when estimating the parameters 7^^, and only the even moments (n + m =even) when estimating the shear 7^, 
since all coefficients with n or m as odd in the background image will have an expectation value of zero. Since deviations 
from zero suggest a second order signal, we are not required to subtract two large modes (e.g. /20 — /02), and thus add 
considerable noise, in order to extract a signal. We discuss a practical inversion technique in the next section. 

3. PRACTICAL PARAMETER ESTIMATION 

While the above relations appear fairly complicated, in reality, inversion is actually quite straightforward. Below, we 
describe a practical technique of inverting a set of shapelet coefficients to yield an estimate of the shear and its derivatives. 

3.1. Measurement and Signal Noise 

As Refregier (2003) notes, if the pixel noise is independent (e.g. Poisson), then the covariance matrix of shapelet 
coefficients will be be simply related by: 



K 



— ""Af^nii 



(33) 



where ctat is simply the Poisson noise of the integrated signal. Though each mode has equal noise, the signal strength 
from the ffist and second order lensing signals are quite different, and it would be helpful to consider in which regime 
each dominates. 

The approximate strength of the second order lensing signal may be inferred from equation (12). Note that all of the 
terms in Dijk are proportional to dj/dr. If the lens is an isothermal sphere (or any other power law relation), then 
Dijk oc "f/Riens, where Riens is the angular distance from the source to lens. Likewise, the integral of XjXkd/dxi over the 
source image will produce a term proportional to the angular scale of the image, itself. We give the "size" of a galaxy 
image as its semi-major axis, Ogai- Thus, the 2nd order lensing signal will have a strength which scales like: 



Ojnm=odd « 7— 



(34) 



which means that unless the source is relatively close to the lens in the image plane, measurement errors can swamp the 
signal. However, it is anticipated (and below, shown), that the intrinsic variance in this signal is significantly smaller than 
the variance in the shear. Moreover, in the particular regime of interest - nearby galaxy pairs, as found in galaxy-galaxy 
fields, this signal will be particularly well-pronounced. 

It should also be noted that, in general, the observed galaxy image is convolved with a PSF. Refregier & Bacon (2003) 
discuss a PSF inversion technique. We do not apply the inversion in the present work, restricting ourselves to large 
galaxies where the impact of the PSF is small. However, we will explicitly perform a PSF inversion in a forthcoming 
paper. It should further be noted that inversion of the PSF gives rise to a non-diagonal covariance matrix, and we thus 
use the more general form of the covariance matrix in our parameter estimation below. 

3.2. Inversion of the Shapelet Coefficients 

In order to compute the lensing coefficients, we need to invert equation (25). In practice, this is somewhat simpler than 
it might initially appear, since, for example, the convergence, k cannot be uniquely estimated from a given shear field (a 
constant k value can be arbitrarily added to the reconstruction). This is known as the "mass-sheet degeneracy," (Kaiser 
& Squires 1993), and, excluding this effect, we may write down a goodness of fit relation as: 
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(7^5",'"' +7*j'S'y 
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7iiinin2m2 



{P-n 



Jn 



:(i) 



(2) 



{l^Si'+J^,jS^^ 



)!< 



(35) 



where the lensing operators implicitly bring power from the primed to unprimcd indices. We have defined, /i„m as the 
"expected" source signal. For n -\- TO=even, this is simply the average of the measured signal (since the universe has no 
preferred direction), and for n + TO=odd, we set this to zero. Likewise, f nm is the best estimate unlensed signal. This is 
subtly different as, though the odd terms are still set to zero, the even terms are the observed coefficients. Since lensing 
(especially in galaxy-galaxy fields) is expected to be small, the best estimate of an intrinsic moment is the measured 
moment itself. 
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Because we have defined the source terms for the odd moments as zero, we may break the above expression into two 
separate terms: 
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^ij^ii h 



(36) 
(37) 



Since the first order term only lenses even moments to even moments, and the second order term only lenses even moments 
to odd moments, the two can be computed independently. In other words, the second-order effects represent an entirely 
new calculation, while the first order x^ can be minimized exactly as described by Refregier & Bacon (2003). 

Once the gradient of the complex shear has been estimated via minimization of ^ , we may relate this to the gradient 
of the convergence using equation (7). We refer to the bending of the observed image as the estimated "Flexion" : 



T = (71,1 + 72,2)1 + (72,1 - 7i,2)j , 



(38) 



where 7^^- is the estimated inversion of the shear derivatives from the x^ minimization. Note that, in principle, T could 
be used to reconstruct a convergence field. In the present worked example, however, noise plays too much of a role. 

3.3. Noise 
With no atmospheric or instrumental effects, the gradient estimator can be related to the true signal via the relation: 



T ^Vk,± (Tyr p ± ajrg 



(39) 



We have already discussed the Poisson measurement noise (crj^,p) above. However, much like in conventional lensing, we 
must also take into account the intrinsic scatter in shapes amongst real galaxies. 

To measure the intrinsic scatter, we observed two complementary regimes. First, we measured the shear and flexion 
variance in the DLS fields themselves. Secondly, we examined the shear and fiexion in two clusters from the HST archive, 
AbcU 665 and Abell 2390, and selected elliptical and spiral galaxies. We found 75 ellipticals and 53 spirals in our claster 
sample which were large and bright enough to classify by eye. No classification was done on DLS galaxies. 

For these samples, we measure the distribution function of ellipticity (shear) and Flexion; the results for clusters are 
shown in Fig. 1 . We define the Flexion in units of the inverse of the semi- major of the observed galaxy. In this way, it is a 
dimensionless, and distance-independent measure of the shape. We find that for both ellipticals and spirals, the standard 
deviation in shear is relatively low, u-y^E = 0.10, and a^,s — 0.13, respectively. 

The scatter in the shear is somewhat lower than that which is normally measured. For example, in the COMBO-17 
sample discussed in Brown et al. (2003), the standard deviation of the ellipticity is approximately 0.25, approximately 
twice as large, and consistent with the estimate from Brainerd et al. (1996). However, taking the same sample, and using 
only those galaxies with r < 17, a^ = 0.16. Since traditional approaches use the KSB technique to invert the PSF for 
small galaxies (an intrinsically noisy inversion), additional scatter is produced in the shear. Since we are interested in 
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Fig. 1. — The fraction of observed galaxies with a particular measured ellipticity (solid) and flexion, ag^i x T (dashed) from HST clusters 
A665, and A2390. Since lensing effects are generally small, and lensing orientation angles are randomized, we suppose this distribution 
to represent the intrinsic galaxy shape distribution of the background sample. Note that while the shear distribution for both spirals and 
ellipticals are similar, ellipticals produce a far smaller scatter in Flexion than do spirals. 
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the intrinsic scatter for the present work, we use our fower estimate of a^. We caution, however, that, even compared to 
bright galaxy sample in COMBO-17, this may be a modest underestimate of the true value. 

For the Flexion, we find a much larger scatter among spirals than ellipticals, with a standard deviation of gf,s = 0.041, 
and op^E — 0.012, respectively. This gives the expected result that early-type galaxies are significantly more regular than 
late types. 

Since sources with lower resolution will be harder to classify, we would like to take a typical Flexion and shear for field 
galaxies. Making the approximation that 60% of a randomly selected background population will be spirals (e.g. Postman 
& Geller 1984), we find a typical shear of a-, = 0.12, and ap = 0.029. These results have a similar to that found in the 
well resolved (but morphologically unclassified) sample from the DLS, yielding, (j^^dls = 0.14, and (Jf,dls = 0.04. 

Using the HST estimates, we may perform a simple estimate of signal to noise. For the first order signal, we find: 
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And, likewise, for the flexion signal. 



since for an isothermal sphere 



since 7 oc 1/R. 

The signals from first and second order effects will be comparable if: 

Rlens ^ -i-icigal (43) 

Thus, at modest separations, the signal to noise from the second order signal will be comparable to the traditional shear 
signal. This ratio will rise to a factor of 10 if the background sample consists entirely of ellipticals. Moreover, it provides 
an entirely new source of information, as it is a direct estimator of the underlying surface density of the lens. 

3.4. Parameter Estimation with a Nearly Circular PSF 

In practice, we do not generally directly estimate cither the shear or the flexion directly from averages over the elliptic- 
ities. Even a perfectly circular PSF will alter the value, but not the orientation, of the estimated parameters. As a result, 
we will generally wish to consider only the relative orientation of the shear or the flexion with respect to the candidate 
lens. 

In practice, this means measuring the Flexion and shear of a large number of background sources which separated from 
their respective lenses by a fixed range of distances. The distribution function of the relative orientation angles of both 
Flexion and shear and then computed. In the absence of lensing, wc would assume these distributions would be drawn 
from a uniform prior. However, lensing will tend to orient the shear perpendicular to the lens, and the Flexion toward 
the lens. 

We must thus relate the distribution function of relative orientation angles, P0, to the induced shear (in first order) or 
flexion (in second order) of the lens. Traditional galaxy-galaxy lensing inversion techniques (Braincrd et al. 1996) yield a 
relation: 

P« = -[l-(7)cos20i(e-i)] (44) 

^ n 

over the domain (fii = [0,7r/2]. The normalization term, (e^^) is exceedingly noisy, and can be unstable for small values 
of ellipticity. Thus, using the parametric model of Braincrd et al. (1996): 

P(e) ex ecxp— A x e , (45) 

Using this relation, we find (e^^ — 14.4) from the a^ discussed above. Note, again, that this is somewhat larger than the 
value of 8 found by Braincrd et al. (1996). 

A virtually identical derivation will yield a relation between the orientation of the flexion: 



p(2) _ 1 



l~{A£-)cOs4,2{{agalJ'y^) 



(46) 



Assuming a similar shape to the shear and Flexion distribution function, and using the ratio a^/ap = 4.1 found above, 
our estimate of the flexion deviation is, {{aiens^)~ ) — 59, over 02 = [0,7r]. 
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4. PROOF OF CONCEPT: THE DLS 

As a proof of concept of this technique, we apply it to the Deep Lens Survey (DLS). The DLS (Wittman et al. 2002) is 
an ongoing deep optical survey of 7x 4 square degree fields, taken on the NOAO 4m Blanco and Mayall telescopes. Each 
field will have an integrated 18ks exposure in the R band, and 12 ks each in the B, V, and z' bands. For the present test, 
we use only R band photometry. While the DLS is designed for measurements of lensing from Large Scale Structure, it 
is also ideally suited for our purposes, as the fields are selected around otherwise empty regions in the sky, and thus the 
primary lensing potential will come from galaxies or large-scale structure only. 

We looked at all 17 of the DLS subfields from the first three public data releases, amounting to a total area of 
approximately 5.4 square degrees. We select as potential lensed sources only those galaxies which have 21.5 < r < 23, 
and which have semi-major axes A > 0.9", since objects significantly smaller than that are not well-resolved into high 
order shapelets. This is a far more conservative cut than on a shallower sample than that used by Brainerd et al. (1996), 
however, we wish to re-emphasize that our goal in this work is to deniontrate the detectability of the Flexion signal. 
Thus, taking as potential lenses those galaxies with 18 < r < 21.5. In total, we found 4833 potential pairs (similar to the 
number found in the Brainerd et al. 1996 sample) with separations less than 60" over 17 subfields. 

We then decompose all background galaxies into shapelet coefficients as described above, and estimate the best fit 
Flexion and shear by minimization x^ as in equation (37). Since seeing produces a significant reduction in the magnitude 
of the shear and Flexion, we use only the orientation angles. 

The distribution of relative orientation angles for pairs with separation less than 16" and greater than 5" is shown in 
Fig. 2, along with a best fit curve representing a shear of 7 = 0.006 ± 0.006 and T = 0.0027 ± 0.0012(")-i. We use 
a KS test in all cases to determine best fit parameters. Since the shear and Flexion should be preferentially aligned 




Fig. 2. — The distribution of relative orientation angles of the shear (top panel) and flexion (bottom panel) with respect to a potential lens 
within a separation of 16". The lines indicate a best fit to equations (44) and (46), of 7 = 0.006 ± 0.006 and dn/dr = 0.0027 ± 0.0012. 



either perpendicular (shear), or toward (Flexion) the source, a rotation of 45 degrees for shear and 90 degrees for Flexion 
(known, hereafter, as the "B-field signal") should produce an entirely random signal. Applying these rotations to each 
observed source and computing the corresponding shear and Flexion produces a means of normalizing the error bars for 
our sample, and checking for consistency. 

For example, if we divide a given sample into iVf,i„ bins of lens-source separation, then we expect that the B-field terms 
will yield a signal with a x^ = Nbin — 1, with the errorbars in each bin being proportional to ^/rNi, the number of galaxies 
within that bin. Moreover, since the errorbars associated with the B-field should be the same as for the E-field, we may 
use this test as a means of normalising the E-field errorbars. We linearly correct extend the B error bars (initially assumed 
to be Poisson noise) such that the reduced x^ is 1. It should be noted that this produces approximately a 10% correction 
in both the Flexion and Shear errorbars from those based on an analytic estimate of uniformly sampled orientation angles. 

This appears to be assuming that there is no systematic error contribution to B; however, we are simply applying a 
useful fiction to obtain error bars for E that contain the systematic error contribution. The method does not explicity 
find the level of this systematic error in B, but accounts for the systematic and statistical errors when quoting an E error. 

We then determine an average radial profile for both terms. We plot the average cumulative shear (measured within a 
disk from separations of 5" to r) for both terms in Fig. 3. 

The radial profiles represent cumulative estimates of the shear and fiexion, and thus the errorbars are strongly coupled 
within a given plot. At highest significance (around 10" for the shear, and 20" for the flexion), the shear and flexion 
result in a 1 and 2ct detection, respectively. The shear estimates are also consistent with estimates from other researchers, 
notably Brainerd et al. (1996) who flnd a shear of 0.0055 for separations of Runs < 20". 
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Fig. 3. — The cumulative estimate of the shear (top panel) and the flexion (bottom panel), of the first three releases from the DLS. Each 
solid point represents the integrated (from Rgai = 5") E-field shear or flexion as described in the text. Note that the errorbars on each are 
strongly correlated. The figure is plotted such that positive shear represents source galaxies aligned tangentially to the lens, and positive 
values of the flexion actually represent decreasing surface density as a function of radius. The squares represent the B-field (rotated) signal, 
and are consistent with no signal. This cross-signal measurement demonstrates that our observed signal could not arise randomly. 



In order to fit the data to physical parameters, we select a simple isothermal sphere model for the lens, and assume 
that all lenses are drawn from the same population. We then expect the relation (Bartelmann & Schneider 2001): 
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(47) 



(48) 



where Dig is the angular diameter distance from source to the lens. We set this ratio to be 0.5 for our discussion since 
lensing signals are typically maximized around this ratio (e.g. Bartelmann & Schneider 2003), and will normalize the 
results to this value. Note that an error in this ratio is systematic, as it effects both our model from the flexion and from 
the shear identically. 

We fit the isothermal sphere model in Fig. 4. Note that we may fit a velocity for either the shear or flexion estimate. 
For the shear, we find a fit of {v^ = I0ltl\km/ s){Ds/ Dis/Q.bf-^ , and for the flexion v^ = {209tllkm/s)(Ds/Dis/0.5f-'^. 
These two fits may be combined to provide a best fit of (201 ± llk'm/s){Ds/Dis/0.5)'^'^. While here we are quoting only 
the random error, there may also be a systematic error based on errors in the assumed distribution function of cllipticities 
and flexions. The shear estimate produces a significantly lower velocity than that found by Brainerd et al. (1996; 220 
km/s), though the two estimates straddle the lower velocity estimate of 135 km/s found by Hoekstra et al. (2004) in the 
Red Cluster Survey. Nevertheless, it should be noted that the shear and Flexion estimate within our own sample are not 
statistically consistent. There are several possibilities. 

First, an isothermal sphere may not be the best model. Consider a circularly symmetric power-law density field, 



K{r) = Ar^'^ , 
where 77 = 1 for an isothermal sphere. The Flexion for such a lens will be: 

Likewise, symmetry dictates that the shear is: 

I7I = ^r) - K{r) = K(r) 

2 — 1] 



Thus: 
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Fig. 4. — As in the previous figure, we plot the estimates of the shear and flexion of galaxies from the DLS. However, here, we plot the 
estimate of the differential shear from annuli of 5". We then fit the curves to an isothermal sphere, as described in the text. Using shear, we 
find Vc = 109 km/s,, while using Flexion, we find Vc = 221 km/s. Again, the squares show the results of rotating the sources by 45 degrees 
for shear and 90 degrees for Flexion. 

Thus, Flexion which indicates a very high velocity dispersion compared to the shear may, in fact, be indicative of an 77 < 1. 
In other words, such a discrepancy may be caused by a density distribution dropping off more slowly than isothermal. 

Secondly, and more likely, recall that though the statistical significance of the Flexion signal can be determined internally 
from the data, the normalization of the signal must be estimated from the intrinsic distribution, in the case of a KS test. 
However, if we've overestimated the intrinsic variance of the Flexion then we will overestimate the Flexion as well. Consider 
that for fixed distribution of orientation angles, the normalizations of the curves may be related as: 



"c.Flexion (^F 
— 7; CX 



V 



(53) 



c. shear 



Thus, the two estimates may be reconciled if the ratio of apja^ is reduced by a factor of (107/221)^ — 0.23. If this is 
the sole source of the discrepancy, then the Flexion may dominate the signal over shear on larger scales than originally 
suggested by equation 43, yielding the new relation: 

^lens ^equality — ^ * -^^gal V^^/ 

as the crossover scale between Flexion dominance and shear dominance. 

5. FUTURE PROSPECTS 

We believe that we have convincingly laid out and demonstrated the feasibility of interpreting new information from 
weak lensing fields by extending that analysis to second order. However, the present work leads to a number of additional 
areas of investigation both theoretical and observational. 

First, the present work does not incorporate inversion of the Point Spread Function. It is therefore clear that we are 
presently unable to extract a significant signal from galaxies which are on order the same size as the PSF. As a result, 
we are forced to remove many potentially lensed galaxies from our sample despite the fact that information could be 
potentially extracted from them. We will address this issue in future work. 

Secondly, the present work has actually thrown away some information with regards to measurements of the shear 
derivatives. Since we have four derivatives, and reduce that data to a 2-vector (J-"), we have, in essence, thrown away 
information. In principle, those derivatives could yield, in addition to dn/dr, a measurement of the variations of, say, the 
shear as a function of radius. Since for an isothermal sphere, these numbers are identical (up to a minus sign), this could 
potentially be an important test for isothermality. 

As a complimentary effort, we will apply this technique to space-based data. Of particular use are the Medium Deep 
Survey (Ratnatunga, Griffiths, & Ostrander 1999) and GOODS (Giavalisco et al. 2004), which will provide the opportunity 
to measure more precisely the intrinsic distribution of flexions in background galaxies. Moreover, since the PSF is much 
smaller, we will be able to directly measure the component of the shear and flexion parallel to the lens-source displacement. 

In summary, we have detected the second-order lensing eff'ect, "flexion" . This effect will be of great value on its own 
and in conjunction with first-order lensing, in order to measure the properties of galaxy halos; it affords direct, local 
information on the gradient of the halo density. 
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